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ABSTRACT 

The planetary nebula populations of relatively nearby galaxies can be easily observed 
and provide both a distance estimate and a tool with which dynamical information 
can be obtained. Usually the requisite radial velocities are obtained by multi-object 
spectroscopy once the planetary nebulae have been located by direct imaging. Here 
we report on a technique for measuring planetary nebula kinematics using the double- 
beam ISIS spectrograph at the William Herschel Telescope in a novel slitless mode, 
which enables the detection and radial velocity measurements to be combined into 
a single step. The results on our first target, the Sab galaxy NGC 4736, allow the 
velocity dispersion of the stellar population in a disk galaxy to be traced out to four 
scale lengths for the first time and are consistent with a simple isothermal sheet model. 

Key words: galaxies: individual: M94, N4736 - galaxies: kinematics and dynamics 
- instrumentation: spectrographs - planetary nebulae: general - techniques: radial 
velocities 



1 INTRODUCTION 



The outer kinematics of galaxies have played a crucial role in 
our understanding of their structure. The dark matter halos 
are most important there, so that conclusions about their 
shape, mass and extent may be drawn that are less depen- 
dent on assumed mass-to-light ratios of the observed stars. 
Most of the angular momentum resides at large radii, and 
relaxation times are longest there, possibly enabling echos 
of the formation process to be observed directly. However, 
the required observations are rather difficult. The integrated 
stellar light of a galaxy rapidly becomes too weak at large 
radii to do spectroscopy. In the case of elliptical galaxies, the 
old stellar populations have now in a few cases been probed 
as far as two effective radii (e.g., Carollo et al. 1995; Gerhard 
et al. 1997). 

Some tracers, such as globular clusters and H I emission, 
can be observed at larger radii, but neither provides a reli- 
able tracer of the kinematics of the relaxed, old stellar pop- 
ulation. Moreover systems like SOs and ellipticals generally 
lack an extensive gaseous disk. Fortunately an alternative 
tracer of the kinematics out to large radii has been identified 
by Hui et al. (1993), who showed that the radial velocities 
of a galaxy's planetary nebula (PN) population constitute a 
suitable diagnostic. Planetary nebulae (PNe) appear in the 
post-main sequence phase of stars in the range 0.8 - 8 M@. 
Fortunately, in all but the very youngest of systems the PN 
population is strongly correlated with the older, and there- 



fore dynamically relaxed, population of low-mass stars. This 
statement is true not only because of the statistics of stellar 
formation and evolution, but also because the PN lifetime 
is itself a strongly decreasing function of progenitor mass 
(Vassiliadis & Wood 1994). 

PNe emit almost all of their light in a few bright emis- 
sion lines, particularly the [O ill] line at 5007A. There is 
evidence that the PN [O ill] luminosity function is essen- 
tially constant with galaxy type and metallicity (Jacoby et 
al 1992), so that the observed bright-end cutoff magnitude 
M* (Ciardullo et al 1989) represents a 'standard candle' 
with which distances can be determined. At a distance of 
10 Mpc, this cutoff corresponds to a flux of 2.0 x 10~ 16 
ergcm _2 s _1 , making PNe within one dex of this limit easily 
detectable in one night with a 4-m telescope. As a rule of 
thumb, approximately 100 such PNe are found to be present 
per 10 9 Lq of B-band luminosity (Hui 1993) , so they are seen 
in sufficient number to study the kinematics of the stellar 
population of the host galaxy. 

The usual approach in making such kinematic stud- 
ies has been to identify the PN population by narrow-band 
imaging and then to re-observe the detected PNe spectro- 
scopically to obtain radial velocities. However, other strate- 
gies exist that avoid the need for several observing runs. 
For example, Tremblay et al. (1995) used Fabry- Perot mea- 
surements in a PN kinematics study of the SB0 galaxy 
NGC 3384. In this paper, we describe a novel alternative, 
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narrow-band Hccfilter narrow-band 0[III] filter 




undispersed Ha image dispersed 0[III] image 



Figure 1. Diagram showing schematically how one can use a 
dual-beam spectrograph to study PN kinematics. The narrow- 
band images in the two arms allow one to identify the PNe and 
measure their line intensities in both the [O III] 5007A and Ha 
lines. The dispersive element in the blue arm shifts each PN image 
by an amount proportional to its redshift. In the red arm the 
grating has been replaced by a mirror. 



based on slitless spectroscopy, and discuss its application to 
the Sab galaxy M94. 



2 DETECTION AND KINEMATICS OF PNE 
THROUGH SLITLESS SPECTROSCOPY 

Our method for obtaining the kinematics of PNe is out- 
lined in Figure [[]. The galaxy under study is imaged through 
narrow-band filters around the two strongest emission lines 
in a typical PN spectrum, Ha and [O ill]. The Ha image is 
recorded directly, but the [O ill] light is dispersed. Compari- 
son of the dispersed and undispersed images then allows the 
kinematics of the PNe to be measured, without prior knowl- 
edge of the location of the PNe. Two modes of analysis are 
possible: 

Dispersed/Undispersed Imaging (DUI): In the dis- 
persed blue arm the PNe will be visible through their [O ill] 
emission as point sources, displaced from their 'true' posi- 
tions by an amount related to their radial velocity, against 
a background of dispersed galactic light. The red arm will 
detect the PNe through their Ha emission, along with any 
other objects with line or continuum emission in the pass 
band of the filter. Assuming that the PNe can be unam- 
biguously identified in the Ha image, their position in the 
[O ill] image will give the radial velocity. We chose to dis- 



perse the blue rather than the red light since PN have a 
higher flux at [O ill] than at Ha, and gratings are less effi- 
cient than mirrors. 

Counter-Dispersed Imaging (CDI): The method of 
counter-dispersed imaging was described in an earlier pa- 
per (Douglas & Taylor 1999). In this mode, pairs of dis- 
persed [O ill] images are made with the entire spectrograph 
rotated by 180 degrees between exposures. The difference 
in the position of a given PN in the two dispersed images 
again reflects its radial velocity. In this case the undispersed 
Ha image can be used as a consistency check on the derived 
positions of the PNe. 

We have implemented our method on the ISIS medium- 
dispersion spectrograph at the Cassegrain (f/10.94) focus 
of the 4.2m William Herschel Telescope. The slit unit was 
removed during the observations, and the [O ill] and Ha 
light paths were separated with a dichroic before being 
passed through appropriate narrow-band filters. The filters 
(A5026/47 and A6581/50) were custom-made for this project 
in order to exploit the full 4x1 arcmin field of the instru- 
ment in slitless mode, and to give adequate velocity cover- 
age. We used a 1200 g/mm (first order) grating in the blue 
arm. Both arms contained 1024 2 -pixel Tek CCD detectors. 
Thus, we obtained dispersed images in the blue (calibration 
showed the dispersion to be about 24 kms -1 per pixel) and 
simultaneous direct images in the red. DUI observations are 
accomplished in a single exposure; CDI requires two expo- 
sures. 

The overall efficiency of this setup, including telescope, 
instrument, filter and CCD, was found from observations of 
a standard star (Feige 34) to be 14% in the blue and 20% 
in the red for air mass of 0. We therefore expected to detect 
~ 2.7 [O in] photons per second from the brightest PN when 
viewed at 6 Mpc. Dark sky (V=21.4) would produce ~2 
counts per arcsec 2 per second, and the background light of 
the galaxy 1-5 counts per arcsec 2 . A reasonable goal is to 
obtain a 4<t detection of the PN population over the top 
decade (2.5 mag) of the luminosity function. This requires ~ 
4 hours of integration if the seeing conditions are of the order 
of one arcsec. The required integration time is approximately 
proportional to the square of the seeing. 



3 OBSERVATIONS 

The observations were carried out on 1997 April 11 and 12. 
For this pilot project NGC 4736 (M94) was chosen. It has a 
large angular size and is at a distance of 6 Mpc (Bosma et al 
1977). The galaxy shows some peculiar morphological fea- 
tures, most notably an inner and an outer optical ring with 
radii of 1 and 5.5 arcmin respectively. The stellar light is too 
faint for direct optical spectroscopy outside 1 arcmin, and 
our goal was to measure PN kinematics to three times this 
radius. Key parameters of NGC 4736 are listed in Table |l|, 
and the observing log is given in Table ^. 

Two fields were observed, 3 arcmin west of centre on 
the major axis, and 3 arcmin north of centre on the mi- 
nor axis. The major axis was observed at two orientations 
(allowing CDI mode) while the minor axis was observed in 
one orientation only. We took the major axis position an- 
gle to be 90°, as would seem appropriate from the relevant 
isophotes (see Figure^). Total integration times were 6.0hrs 
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Table 1. Galaxy parameters 



Name 


NGC 4736 (M94) 


Position (J2000) 


12h50m53.061s 




+41d07ml3.65s 


Hubble type 


SA(r)ab 


VHolio 


308 ± 1 km/s 




8.99 


Angular size 


11.2 X 9.1 arcmin 


Inclination 


M 35 ° 


Scale length 


57 i 10 arcscc 


Position angle 


105° 



Table 2. Log of our observations at the WHT, La Palma Obser- 
vatory. PA is the position angle of the spectrograph and T the 
total integration time (where two values are given these refer to 
the red/blue arms respectively). 



Field PA T UT Airmass 











April 11 






HD66637 




90 


60 




20:35 




1.06 


Peige 34 




90 


2/10 




21:15 




1.07 


N4736 - major 


axis 


90 


16,453 


21:20 


- 2:13 


1.35 - 


■ 1.07 


PN 49+88.1 




90 


5 




2:03 




1.06 


N4736 - major 


axis 


270 


4,868 


2:36 


- 4:52 


1.12 - 


■ 1.54 










April 12 






HD66637 




90 


0.1/10 




20:29 




1.06 


Feige 34 




90 


3/60 




20:57 




1.09 


N4736 - minor 


axis 





10,800 


21:02 


- 0:08 


1.42 - 


1.03 



on the Western field (4.6hrs in one orientation, and 1.4hrs 
with the spectrograph rotated by 180 degrees), and 3hrs on 
the Northern field. The observing conditions were close to 
photometic with seeing, as judged from stellar images in the 
red arm, less than 1.1" at all times. We also observed a flux 
standard star for photometry and a Galactic PN as a radial 
velocity reference. A second star was observed as a spectral 
reference. 

The custom-made [O m] and Ha filters had a cen- 
tral wavelength and peak transmission of 5026A/0.823 and 
6581A/0.915, respectively. The nominal FWHM was 47A 
and 50 A, while the effective photometic bandwidth was eval- 
uated graphically and found to be 38. 7A and 45. 7A, respec- 
tively. 



4 DATA REDUCTION 
4.1 Calibration 

The dispersion in the blue arm was measured by inserting a 
slit and using an arc lamp, and found to be 0.3992A/pixel. 
The spectrum of the star HD66637 was then observed 
through the same slit and wavelength-calibrated. Subse- 
quent observations of the same star at numerous positions 
in the field (after removal of the slit) established that the 
dispersion could be taken as constant over the field. The 
combination of these observations with the undispersed red 
arm positions of HD66637 gave an unambiguous solution for 
the transformation between objects in the red (direct image) 
and blue (dispersed image). (Note that with this technique 
the radial velocity of the star does not enter into the calcu- 
lation.) 



To check the zero point of the velocity scale, we moved 
the telescope from the reference star to the Galactic plane- 
tary nebula PN 49.3+88.1, for which the heliocentric radial 
velocity is listed as —141 kms -1 (Schneider et. al 1983). In 
ten pointings over the field of the spectrograph we measured 
— 136.0 ± 3.2 kms -1 , in agreement with the calculated ob- 
servatory frame redshift of —133 kms" 1 . Unfortunately we 
discovered later that at certain telescope orientations the 
flexure is large enough to cause significantly larger errors. 
However, such flexure only introduces an offset in absolute 
velocity, and does not compromise our ability to study a 
galaxy's internal kinematics. In order to derive an absolute 
calibration for the velocity scale, we later obtained a long 
slit observation of two of the objects detected in this analysis 
(see § |). 

The spectrograph field of view with the slit unit re- 
moved consists of an approximately unvignetted area of 
about 4'x 1', but we obtained useful data outside this re- 
gion. Correcting the observed fluxes for the vignetting is only 
straightforward for the red (undispersed) arm. In the blue 
arm the correction is complicated by the fact that the im- 
age is dispersed. The sky flat measured in the blue arm was 
found to closely approximate the aperture function in the 
red arm, shifted by a small number of pixels, transformed 
to blue arm coordinates and then convolved with the filter 
profile at the appropriate dispersion. 

We therefore carried out the complete analysis after cor- 
recting only for the pixel-to-pixel variation of the CCD re- 
sponses, determined in the usual way. Once the PNe were 
identified it was possible to determine their wavelength and 
their positions in the aperture prior to being shifted by the 
spectrograph, so that the [O ill] magnitudes could then be 
corrected analytically both for the filter response, which was 
fitted with a polynomial, and for vignetting. 



4.2 Object identification 

Scripts based on iraf procedures were used for all of the 
data reduction. 

Due to spectrograph flexure, individual exposures 
needed to be aligned before being added. This was accom- 
plished by a simple shift (at least one PN was visible in 
each individual exposure). The images were then combined 
by computing a weighted and scaled median. A spatial me- 
dian filter was applied to the combined frame, and the result 
subtracted from the original image to yield a field of unre- 
solved objects against a background with a mean of zero. 

The PNe and any other point sources in the red and 
blue images were extracted by two methods: 

(A) Blinking and hand-tagging, followed by a PSF-fitting 
step to evaluate the shape and size parameters. For DUI 
mode observing it was usually found to be easier to search 
the range of possible (red) coordinates corresponding to an 
object seen in the blue (see Fig. ^). 

(B) an automated procedure based on object lists generated 
with daophot. A 2D-gaussian fit to each detected image was 
used to select PN candidates. The FWHM of a candidate 
was required to be within a small range of that of the seeing 
disk (PNe are unresolved) and the axial ratio close to the 
value 1.28 expected from our instrumental configuration (the 
ellipticity arises from the anamorphic effect of the grating). 



© 1999 RAS, MNRAS 000, §-| 



4 N.G. Douglas, J. Gerssen, K. Kuijken and M.R. Merrifield 




Figure 2. An example of the 'postage stamp' method for detect- 
ing red counterparts. The boxes on the right each contain one 
point source identified in the blue frame. The boxes on the left 
side show the corresponding parts of the red frame. Any red coun- 
terpart must be located somewhere within the box since the size 
of the box is such that it encompasses the velocity range spanned 
by the filter. The horizontal location of the counterpart is a mea- 
sure of the object's radial velocity. In the vertical direction we 
require the positional agreement to be within one pixel. 

The object lists were then correlated to search for potential 
PN image pairs. 

4.3 False detections 

As well as PNe, our observations will also detect other ob- 
jects having a predominantly line-emission spectrum near 
the [O ill] line. In the case of a bright spiral galaxy like M94 
it is obvious that H II regions will mimick PNe. As they be- 
long to a younger population than the PNe their inclusion in 
the analysis would lead to an underestimate of the velocity 
dispersion. Although it is tempting to use the [O ill] /Ha line 
ratio as a discriminant (it tends to be larger in PNe) this 
is not sharp enough to eliminate the H II regions without 
eliminating a fair fraction of the PNe. Unless a better dis- 



criminant can be found, this source of contamination may 
ultimately restrict the viability of our technique to the el- 
liptical and SO galaxies for which it was devised. However 
in the particular case of M94 the problem is manageable 
because: 

(i) M94 is sufficiently close that a large fraction of bright 
H II regions would be spatially extended and thus eliminated 
by one of the object selection criteria; (ii) the H II regions in 
M94 are mostly confined to features associated with spiral 
arms - one of these passes through the eastern extremity 
of our major axis field and all objects there were ignored. 
The rest of the field, and the minor axis field, are devoid 
of recognised H II agglomerations; (iii) the total number of 
unresolved objects detected in [O ill] is close to the number 
expected for PNe on the basis of their luminosity specific 
density as seen in other galaxies. 

None of these factors eliminates the problem entirely, 
but their combined effect is such that we feel that con- 
tamination is negligible. Even if 10% of the PNe had been 
misidentified the effect on the calculated velocity dispersion 
would be at most 5%, well below other sources of error. 

High redshift galaxies, in which the Lya line is shifted 
into the [O ill] passband, form another potential source of 
contamination when long integrations are made of extended 
halos (Freeman et al. 1999). This is also not a significant 
issue in the present case. 

4.4 Comparing Spectral Modes 

We identified PNe in NGC 4736 along the minor axis with 
DUI mode observations and along the major axis using CDI 
(with DUI mode data being redundant). Although CDI re- 
quires two distinct integrations, we found that data obtained 
from CDI results in a higher number of detected PN per 
unit integration time. For visual identification using blink- 
ing, CDI is considerably easier since the two images have 
the same plate scale and similar properties with respect to 
sky noise and confusion. Therefore to illustrate the numeri- 
cal superiority of using two dispersed [O in] images we rely 
on the automated search results for the major axis observa- 
tions only. 36 PNe were identified in this field from matching 
2.5cr sources in the counterdispersed [O ill] images. The lim- 
iting factor here was the shorter integration time with one 
of the two spectrograph orientations: had both integration 
times been equal we would presumably have found yet more 
PNe. By comparison, only 24 PNe were detected from the 
DUI mode analysis of the same field. Therefore, we conclude 
that a significant number of PNe are too faint in Ha for the 
DUI mode to detect them. 



5 LONG-SLIT DATA 

It has been mentioned that instrument flexure, particularly 
between sets of observations in CDI mode, can lead to an 
uncertainty in the absolute velocity scale. To remove this 
uncertainty we attempted to obtain the velocity of at least 
one object in the major axis field via the Willian Herschel 
Telescope service data program. As the PNe are faint, with 
effective V-band magnitude of around 25, the slit had to 
be positioned 'blind' on the basis of the position computed 
from the dispersed data. Fortunately the astrometry does 
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Figure 3. Digitised Sky Survey image of NGC 4736 (North at 
top) with the observed fields marked with boxes and the PNe 
shown as open circles. 



not depend on velocity, but only on correct identification of 
PN pairs, and on the centroiding of the dispersed images of 
stars in the field. 

We requested a spectrum using ISIS in long-slit mode 
and with the slit at a position and PA chosen such as to 
fall across two objects. This provided a good test of the 
astrometry. The service observation was attempted on 1999 
July 28. Both objects were acquired, and their separation 
along the slit agreed with that calculated. One of the objects, 
suspected of being an H n region, was confirmed as such. 
The radial velocities obtained had an internal error of about 
10 kms -1 , as judged from the values obtained from different 
lines, and were used to calibrate the major axis data (Table 



6 RESULTS 

The PNe identified along the major axis are listed in Ta- 
ble ^ and those along the minor axis in Table ^ Their po- 
sitions are also shown in Fig |[ As suggested by the suc- 
cessful long-slit experiment, the positional uncertainty is of 
the order of 1-2 arcsec. The internal error in the velocities 
is approximately 10 kms -1 . The minor axis velocities have 
an offset that has not been determined, but the values as 
presented have an average velocity near systemic, as would 
be expected. The derivation of velocities from DUI mode 
data, as was used for the minor axis, is in fact much less 
sensitive to flexure. However, in order to reduce the number 
of candidate objects in the red image, only radial velocities 
between 200 and 500 kms -1 were searched for, so this table 
has to be used with caution. 



Table 3. Fifty-three point sources in the major axis field. Rvel 
is the relative (heliocentric) radial velocity measured from the 
slitless spectroscopy plus an offset established by an additional 
calibration (see §H). rri5007 = -2.5 log(F) - 13.74 is the appar- 
ent V-band magnitude of the objects, where F is the flux in 
ergc m - ^ -1 . The magnitudes have been corrected for vignetting 



(see f 4.1). The vignetting correction could not be reliably applied 
for two objects as they were too close to the edge of the aperture, 
and their magnitudes are not given. FWHM (in pixels) refers to 
the major axis of the 2D-gaussian fit. 



da ( T9nnn^ 


Dec (.12000) 


Rvcl 


m 5007 


C vv X11V1 


12:50:49.46 


41 


07:09.5 


45!) 


23.44 


3.52 


12:50:49.42 


41 


07:40.3 


472 


23.18 


3.59 


12:50:49.35 


41 


07:17.5 


463 


23.48 


3.31 


12:50:48.93 


41 


07:51.2 


483 


22.54 


3.29 


12:50:48.62 


41 


07:22.9 


444 


23.18 


3.31 


12:50:48.40 


41 


06:17.3 


309 


24.80 


3.24 


12:50:48.53 


41 


07:45.8 


450 


24.63 


3.54 


12:50:47.97 


41 


07:10.8 


508 


25.46 


3.9 


12:50:47.87 


41 


07:09.8 


378 


25.78 


2.79 


12:50:47.73 


41 


06:33.2 


296 


24.56 


3.9 


12:50:47.61 


41 


07:00.2 


345 


24.45 


3.61 


12:50:47.59 


41 


06:56.9 


391 


25.33 


3.38 


12:50:47.50 


41 


06:26.2 


525 


25.52 


3.5 


12:50:47.17 


41 


06:27.4 


42!) 


25.34 


3.54 


12:50:46.92 


41 


06:05.5 


413 


25.14 


3.45 


12:50:46.93 


41 


06:35.3 


378 


25.41 


3.47 


12:50:46.32 


41 


06:18.3 


357 


25.60 


3.05 


12:50:46.40 


41 


07:46.8 


466 


24.95 


4.38 


12:50:46.13 


41 


06:11.2 


380 


24.24 


3.94 


12:50:46.17 


41 


08:13.6 


491 


24.34 


3.57 


12:50:45.27 


41 


06:36.8 


478 


24.91 


4.37 


12:50:45.11 


41 


06:27.8 


308 


25.30 


3.73 


12:50:44.65 


41 


06:51.3 


369 


24.82 


3.66 


12:50:44.26 


41 


07:39.2 


432 


25.78 


2.67 


12:50:43.79 


41 


06:11.4 


353 


24.95 


3.33 


12:50:43.76 


41 


06:48.3 


348 


25.18 


2.79 


12:50:43.29 


41 


07:54.9 


457 


22.59 


3.52 


12:50:43.03 


41 


06:07.7 


485 


24.91 


3.17 


12:50:42.95 


41 


07:49.5 


473 


26.39 


2.7 


12:50:42.66 


41 


06:11.8 


405 


22.65 


3.92 


12:50:42.58 


41 


06:36.6 


428 


25.00 


4.34 


12:50:42.48 


41 


07:52.8 


4 2 5 


24.88 


3.5 


12:50:42.19 


41 


08:15.0 


470 


24.58 


3.07 


12:50:41.40 


41 


07:01.5 


433 


24.98 


3.5 


12:50:41.17 


41 


06:21.0 


413 


25.21 


3.47 


12:50:41.17 


41 


07:46.4 


339 


24.26 


4.7 


12:50:40.87 


41 


05:57.5 


414 




3.47 


12:50:40.74 


41 


06:11.7 


410 


25.53 


4.37 


12:50:40.55 


41 


07:46.4 


417 


24.72 


3.61 


12:50:40.26 


41 


08:00.7 


400 


25.33 


4.04 


12:50:40.15 


41 


06:51.3 


383 


24.74 


2.77 


12:50:38.62 


41 


07:04.9 


40!) 


25.41 


3.71 


12:50:38.67 


41 


08:05.0 


391 


25.35 


4.2 


12:50:38.55 


41 


06:44.4 


418 


25.00 


2.89 


12:50:38.02 


41 


06:20.9 


423 


25.54 


3.68 


12:50:37.90 


41 


08:06.1 


404 


25.14 


3.05 


12:50:36.34 


41 


06:48.0 


430 


25.65 


3.07 


12:50:35.24 


41 


06:10.8 


425 


24.93 


3.45 


12:50:34.75 


41 


06:18.3 


408 


24.87 


3.31 


12:50:34.15 


41 


05:58.2 


379 




2.91 


12:50:25.18 


41 


06:16.4 


415 


25.50 


4.3 


12:50:48.71 


41 


06:32.7 


378 


25.18 


5.45 


12:50:42.79 


41 


07:48.9 


454 


25.22 


3.54 



6.1 Luminosity function 

We placed a premium on detecting as many PNe as possi- 
ble, even in the partially vignetted region of the instrument. 
Considerable corrections have been applied, and the mag- 
nitudes should therefore only be taken as indicative. The 
luminosity function of the objects detected is presented in 
Fig ^. The bright-end cutoff for the assumed distance of 
6 Mpc (m* = 24.4) is indicated. At the faint end the lumi- 
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Table 4. Fourteen point sources in the minor axis field. The 
coordinates are derived directly from the position in the red image 
(see text), rgp^gi gives the ratio of the fluxes (in ergcm _2 s — 1 ) 
seen in the red and blue images. Other symbols are as in Table jj| 
except that the radial velocities have not been checked against an 
external reference and may have a small common offset. 



RA (J2000) 


Doc (J2000) 


Rvel 


m 5007 


Ha 
TO~raT 


FWHM 


12:50:52.02 


41:11:22.3 


330 


26.26 


1.26 


3.36 


12:50:53.68 


41:10:15.8 


263 


24.37 


0.38 


3.17 


12:50:50.03 


41:09:22.2 


I'M 


26.56 


1.11 


3.71 


12:50:50.16 


41:09:01.5 


311 


25.67 


0.82 


4.04 


12:50:51.57 


41:08:58.3 


270 


24.74 


0.62 


3.1 


12:50:56.25 


41:08:42.6 


265 


25.54 


0.65 


4.74 


12:50:52.22 


41:08:35.6 


268 


24.35 


0.48 


3.54 


12:50:48.81 


41:08:25.3 


367 


26.21 


1.57 


2.72 


12:50:57.74 


41:08:19.7 


252 


25.35 


0.65 


4.25 


12:50:46.14 


41:08:13.2 


408 






3.61 


12:50:49.62 


41:08:04.5 


363 


24.84 


1.28 


5.02 


12:50:52.93 


41:08:02.5 


496 


25.63 


2.47 


4.46 


12:50:53.86 


41:07:56.1 


310 


23.35 


2.61 


3.31 


12:50:51.57 


41:07:52.1 


367 


23.62 


8.87 


3.76 



nosity function is, of course, significantly incomplete while 
at the bright end some objects are brighter than the cutoff. 
The latter are probably H II regions and have therefore not 
been included in the analysis of the kinematics. 

The number of PNe found in the major axis field 
(CDI mode) is in rough agreement with predictions. From 
the basic data on NGC 4736 compiled by Mulder (1995) 
(m B = 8.58, D = 6.0 Mpc) we have L B = 2.07 x 10 10 L Q , 
and on the basis of the results of Hui et al (1993) the ex- 
pected number of PN in the top decade of the PNLF in 
M94 would therefore be around 2000. Mulder also found the 
galaxy to be fairly well-fitted by an exponential disk with 
scale length h =57". The region we examined includes 0.029 
of the light of such an exponential, which should therefore 
include 59 PNe in the brightest decade. This number com- 
pares well with the 53 actually found, though the agreement 
may be somewhat fortuitous given the incompleteness at the 
faint end. 



6.2 Rotation curve 

In Fig ^ the line-of-sight velocities of the 53 objects in the 
major axis field are plotted as function of radius, after sub- 
traction of the systemic velocity (Table |l|). Flat rotation is 
seen out to the last measured point at almost 5 scale lengths. 
For comparison we overplot the H I / CO rotation curve of So- 
fue (1997), projected into the plane of the sky. For objects 
near a distance of 1 arcmin along the major axis the mean 
velocity is 98 kms -1 , in agreement with the projected gas 
rotation velocity of 103 kms -1 at that point. The uncor- 
rected minor axis data yield a mean velocity of 329 kms^ 1 , 
consistent with the systemic velocity. 



24 

m5007 



Figure 4. The luminosity function of the 64 objects for which 
magnitudes could be determined. The cutoff peak for PNe is 
clearly defined. The objects brighter than this cutoff are prob- 
ably bright H II regions 



K 



a " B 
- —a 



2 4 6 

Distance along major axis (arcmin) 

Figure 5. Line-of-sight velocities of the 53 objects in the major 
axis field plotted as a function of distance along the axis. The 
dotted line shows the gas rotation curve, from Sofue (1997) 



6.3 Velocity dispersion 

Generally, the vertical structure in disks of spiral galaxies 
is reasonably well described by an isothermal sheet approxi- 
mation (van der Kruit & Searle 1982; Bottema 1993). In this 
model, the vertical velocity dispersion is found to follow an 
exponential decline with radius with scale length twice that 



of the surface density. With the additional assumption that 
the dispersion ellipsoid has constant axis ratios throughout 
the disk, one finds that the line-of-sight velocity dispersion 
follows the same decline, independent of the galaxy's incli- 
nation. We tested this using the major axis data. Figure 
shows the velocity dispersion in bins of distance along the 
major axis. Seven objects were eliminated from the kine- 
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matic analysis as their brightnesses suggested that they are 
H II regions. The curve is the least-squares exponential fit, 
yielding central velocity dispersion 111 kms -1 and scale 
length h a = 130 arcsec. These are close to the published 
value for the central stellar velocity dispersion (120 ± 15 
kms -1 ) obtained from absorption- line spectra (Mulder and 
van Driel, 1993) and to twice the photometric scale length 
(2ft = 114 arcsec), suggesting that the isothermal sheet ap- 
proximation is reasonable. 



6.4 Combined Kinematic Model 

The binning in Fig ^ effectively assumes that the PNe all lie 
close to the major axis. In fact, they are located up to one 
arcminute from the axis, at azimuths up to 45° , so a more so- 
phisticated approach is required. We therefore projected the 
PNe on to a thin disk of fixed inclination (35°), giving r, <f> 
coordinates. The nebulae's line-of-sight velocities can then 
be compared with a model consisting of a three-dimensional 
isothermal sheet with a flat rotation curve. This model has 
five parameters, namely the three components of the cen- 
tral velocity dispersion o z , o^, or, the scale length h a , and 
the rotation amplitude. A maximum likelihood method was 
then used to fit the model. We added twelve PNe from the 
minor axis field (Table ^) to help constrain the fit (two were 
excluded from the fit as probable H II regions). 

In practice the data were not adequate to constrain all 
five parameters. Using the canonical relationship cr^ = o R /2 
from the epicyclic approximation (Binney & Merrifield 1998, 
eq. 11.18) and allowing <t z /(Tr to vary over the range 0.2 to 
2.0, we found a robust maximum likelihood solution with 
scale length ft CT — 144 ±30 arcsec, central velocity dispersion 
oios = 120 ± 30 kms -1 , and circular rotation speed v c = 
177 ± 11 kms -1 . These results are consistent with the H I 
rotation speed at 1 arcmin radius (180 kms -1 ) and again 
with twice the photometric scale length. 

It was not possible to constrain the shape of the ve- 
locity ellipsoid with these data - such an analysis would 
require a more complete azimuthal coverage of the galaxy. 
However if we assume 0.5 < o z /or < 1.0 then we infer 
75 < or < 110 kms -1 at one photometric scale length, 
consistent with the trend between rotation speed and disk 
velocity dispersion found by Bottema (1993). 

Thus far we have ignored measurement error in the ve- 
locities, which will tend to increase the measured velocity 
dispersion. This turns out to be a small effect: allowing for 
a lcr error of 10 kms -1 , the fitted dispersion becomes ap- 
proximately 3% smaller and the scale length is unchanged. 



120 




Distance along major axis (arcmin) 



Figure 6. The dispersion in radial velocities computed by bin- 
ning the radial velocity measurements, together with the best-fit 
exponential curve 



PNe are not reliably high enough to allow both spectral lines 
to be used. 

Our pilot experiment was performed on the large nearby 
Sab galaxy M94. It has revealed a PN population in the disk 
whose rotation curve remains flat, and whose velocity dis- 
persion declines radially exponentially, consistent with the 
predictions of a simple isothermal sheet model. PNe were de- 
tected out to five exponential scale lengths, well beyond the 
reach of kinematic measurements based on integrated-light 
absorption- line spectroscopy. The number of PNe detected 
was consistent with expectations. 

The present experiment was limited to two fields in 
this large galaxy. Complete coverage of the galaxy should 
yield around 2000 PNe, and would allow a detailed kine- 
matic model to be fitted, including a determination of the 
axis ratio of the velocity ellipsoid following the technique of 
Gerssen et al. (1997). Obtaining such data for a small sam- 
ple of nearby galaxies in just a few nights of 4-m telescope 
time is a practical proposition. 
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7 CONCLUSIONS 

In this paper we have demonstrated how the kinematics of 
the PN population in a galaxy can be measured by slit- 
less spectroscopy through narrow-band filters with a dual- 
beam spectrograph. We compared two possible modes: dis- 
persed/undispersed imaging, in which a dispersed [O ill] im- 
age is compared to an undispersed Ha image; and counter- 
dispersed imaging, in which two [O ill] images, dispersed in 
opposite directions, are analysed. It turns out that the lat- 
ter method is more effective: evidently the Ha fluxes of faint 
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